SEP an 1949 


AERONAUTICAL 
QUARTERLY 


Volume I AUGUST 1949 Part Il 


CONTENTS 


Displacements of a Linear Elastic System 
Under a Given Transient Load 


Skin Friction in the Laminar Boundary 
Layer in Compressible Flow A. D. Young 


The Measurement of Gas Turbine Com- L. J. Richards 
bustion Efficiency by Gas Analysis and J. C. Street 


A Simplified Theory of “Simple Waves ” A. Ghaffari 


LONDON 
ROYAL AERONAUTICAL SOCIETY 
4 HAMILTON PLACE WI 


| ij 

| f AE 

| 

| 10/- 


THE AERONAUTICAL QUARTERLY 


Editorial Board 
W. G. A. PERRING, F.R.Ae.S. (Chairman) 
Dr. H. ROXBEE COX, D.LC., F.R.Ae.S., F.1.Ae.S. 
Professor S. GOLDSTEIN, F.R.S., F.R.Ae.S. 
Sir RICHARD SOUTHWELL, M.A., LL.D., F.R.S., F.R.Ae.S. 


G. H. DOWTY, F.R.Ae.S. (Chairman of the JouRNAL Committee 
of the Royal Aeronautical Society) 


J. SMITH, C.B.E., F.R.Ae.S., A.M.I.Mech.E. (Chairman of the 
Technical Board and Executive Committee of the Society of 
British Aircraft Constructors) 


Editorial Executives 


Captain J. Laurence Pritchard, Hon.F.R.Ae.S., Hon.F.].Ae.S. 
Joan Bradbrooke, A.R.Ae.S. Dr. D. M. A. Leggett, A.F.R.Ae.S. 


Editorial Referees 


Aerodynamics Instrumental and Aircraft Propulsion 
Electrical Equipment 


. G. E. Bamsto 

Mr. I. BowENn 
. G. W. H. Garpner 
. C. S. Hupson 


tructures 
A. C. Brown 


Fluid Motion 


Mr. S. 
. . R. H. SANDIFER 


RICHARD SOUTHWELL 


els and Oils 


REY 


BERR 


ircraft Accessories 


> 


Radio and Radar 
{c) Stability and Control _ J. F. ATHERTON nae 
L. W. . B. G. Gates 
Mr. S. B. . R. L. Smirn Rose 
. A. M. 
Rosert WaTsON WATT 


C. G. 


G 
. 

M 


Aviation Medicine 
Helicopters and Dr. A. BUCHANAN BARBOUR 


Propellers 
P Meteorology Performance 
Mr. H. Davies 
Mr. D. R. Newman 
Mr. A. H. Yates 


Prof. L. Bamsrow H. L. Cox 
Dr. G. P. Dougias D. J. FARRar 
Prof. A. A. HaLbi H. GRINSTED 
Sir B. MELVILL Jones H. Roxsee Cox 
Mr. D. Kerra-Lucas Mr. A. N. IRens W. S. Hemp 
Dr. N. A. V. Preacy Mr. C. N. Jaques H. B. Howarp 
Mr. W. G. A. PERRING Mr. F. W. MEREDITH E. D. KEEN 
Mr. E. F. RELF Mr. H. C. Pritrcuarp F. A. Kerry 
Mr. R. LEGoETT 
. E. J. Richarps ; . L. Licxiey 
Aircraf Design A. J. SUTTON Pipparp 
E. L. Bass 
J. W. Drinkwater 
F. H. Garner 
W. E. Garner 
Materials 
Prof. L. AITCHISON 
bration and Flutter Dr. N. P. ALLEN 
f. A. R. COLLaR Prof. N. F. Morr 
Mr H. Lee Prof, W. J. DUNCAN Mr. L. ROTHERHAM 
Mr. R. A. FRAZER Dr. H. SUTTON 
Mr. J. HANSON Dr. C. Syxes 
Flight Testing 
Mr. P. A. HuFrTon 
Mr. E. T. Jones 
Mr. S. Scotr HALL 
Mr. H. L. STEVENS 


j 
0 


DISPLACEMENTS OF A LINEAR ELASTIC SYSTEM 
UNDER A GIVEN TRANSIENT LOAD 
by 


D. WILLIAMS, D.Sc., M.I.Mech.E., A.F.R.Ae.S. 
(Structures Department, Royal Aircraft Establishment) 


PART I 


1. INTRODUCTION 


A method is described for obtaining the displacements (and hence the stresses) 
of a linear elastic system under a given transient load in a more accurate and more 
expeditious manner than is possible by the standard method. For simplicity, no 
account is taken of damping forces in describing the method. 


The method is of interest to engineers and others who require to obtain 
numerical answers to problems of the type described above. It depends upon the 
fact that in many practical cases of transient loads the inertia forces in the series of 
normal modes of an elastic system converge rapidly, compared with the series of 
“normal” loadings in terms of which the spatial distribution of the applied load 
can be expressed. 


The first application of the method to a practical problem was made by the 
author in a Royal Aircraft Establishment”) report in which the effects of landing and 
of gust forces on the wing stresses of an aeroplane were investigated. The essence of 
the new approach was subsequently described in a paper read at the Sixth Inter- 
rational Congress of Applied Mechanics (Paris, September 1946), and the form of the 
argument used there is repeated in Sections 2-5. For readers familiar with vibration 
theory the formal presentation of the argument in those Sections requires no 
amplification; but the practical usefulness of the method is so great, particularly for 
aeronautical engineers, that it was thought worthwhile to expand the argument, and 
to show, by means of a simple example, precisely how the saving in labour and gain 
in accuracy is brought about. In the matter of accuracy, it will also be made clear 
why, theoretically, the standard method is incapable of giving an accurate solution 
and why complete accuracy is achieved by the modified approach. 


Paper received December 1948. 


[The Aeronautical Quarterly, Vol. 1, August 1949] 
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2. BASIS OF STANDARD METHOD 
The argument runs as follows : — 


Let a force P (x). F (t) be a transient load applied to a continuous elastic system 
(with therefore an infinite number of degrees of freedom), where P (x) represents 
the spatial distribution (assumed constant with time) of the applied load and F (ft) 
its variation with time. 


It is well known that the function P (x) can be expanded in a series of normal 
loadings each having a distribution identical with _ of the inertia loading in the 
corresponding normal mode. Thus let 


P(x)=L, (x)+L, (x)+..... +L, (x)+...... 
and therefore 


3. STANDARD PROCEDURE 


The standard procedure is to calculate the displacements caused by each of 
the terms on the right hand side of (la) independently of the others, beginning with 
F (t) L, (x) and continuing one by one until it is considered that a close enough 
approximation to the total effect has been reached. 


The disadvantage of this procedure is seen immediately when, for example, 
only the first two or three resonance modes of the system are excited, the remaining 
modes having displacements corresponding only to the static application of the 
corresponding loadings. These remaining modes have to be considered up to a 
term F (t) L, (x) where r is high enough to make 


M- 


L,, (x)=P (x) 


n 


to a reasonable degree of accuracy. 


If P (x) is a concentrated load r is likely to be fairly large, making it necessary 
to evaluate a correspondingly large number of terms. 


4. SUGGESTED PROCEDURE 


To obviate the difficulty indicated above we note that, in the absence of the 
corresponding inertia forces, the displacement of the system in the n' mode at 
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every instant is equal to the statical displacement appropriate to the nm” normal 
loading, i.e. to the force F (t). L, (x) at that instant. If we add to this the displace- 
ment corresponding to the appropriate inertia loading f, (t). L, (x) we shall have the 
resultant displacement of the system in the n mode, the function f(t) representing 
the time variation of the inertia loading L, (x) and being different for each mode. 


The total resultant displacement at every instant is therefore that appropriate 
to the loadings: — 


. 
F (t) La(x)+ fr) La) 
n=l i=l 
or to P (x). F (t)+ fu (0). Ln (2). 


Now in many practical cases the second term in (3) converges rapidly and can 
be evaluated with little difficulty. The first term, being the original applied load, 
is already known and its effect (as a static force) can be added immediately to that 
of the second term to give the complete behaviour of the system. 


5. EVALUATION OF INERTIA FORCES 


The inertia forces are readily obtained by the use of Duhamel’s Integral. In 
the n mode, for example, the Equivalent Static Force is given by 


t 


=p, @) F (t,)sin p, (t—t,)dt, . (4) 


where (E.S.F.), is the static force in the n'* mode that gives the same displacement 
as the actual dynamic force L,, (x). F (t) at every instant t from the beginning of 
force application, and where P,, is the circular frequency in the same mode. 


Now, at every instant the Applied Force (A.F.), (i.e. L, (x) F (t)) is resisted by 
an elastic force equal to (E.S.F.),, but opposite in sign, and an Inertia Force (I.F.), 


125 


F(t) 


D. WILLIAMS 


F(t,) 


ANO 


Fig. 1. 
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(already referred to above as f, (t). L, (x) ), and therefore 


(L.F.), = —(A.F.), + (E.S.F.),, 


t 


=L, | F (t,) sin p, | : : (5) 


This gives, on integration by parts, 


(LF.),=L, —F F (t)- cos p, (t—t,) F (t,)dt, | 


t 
=~ L, (x) | cos pa (t—1,) F (t,) dt, (6) 


where F(t) =F (t,)], and the legitimate assumption is made that the initial 
value of F (t) is zer0. 


If Fig. 1 represents the curve of F(t) and Fig. 2 the curve of F (¢) against 
time f¢ it is clear that as p, increases, the value of the integral in (6) decreases (unless 
sharp changes of slope take place in time intervals comparable with 27/p,). 


In many practical cases (I.F.),, becomes small enough to be neglected after the 
first few terms so that 


(LF.),= = fn (t). Ln (x) = (LF.), . (7) 


(where r is small), is a series that is readily evaluated. 


Resultant effect 

The resultant displacements, and hence the resultant stresses in the system, are 
therefore found at every instant from equation (3) by adding the static displacements 
produced by the inertia loadings to that produced by the actual applied load. 


6. SIMPLE EXAMPLE 


The convenience of the method will be appreciated if it is applied to the easy 
example represented by a uniform beam, simply supported at the ends, and 
subjected to a dynamic load applied at the mid-section and varying harmonically 
with time. Figs. 3(a@) and (b) together define completely the conditions of the 
problem, (a) giving the spatial distribution of the applied load, and (5) its variation 
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with time F (t), which has here the form sin p.f. We shall consider the behaviour 
of the beam from zero time—representing the onset of the load—to time 7/p, when 
the force is removed after passing through its maximum value Q at time z/2p,. We 
shall assume that the fundamental circular frequency p, of the beam is given by 


p.=p/4 . : (8) 


from which, in the present case, the higher frequencies can be written down. As, 
by symmetry, modes that have a node at the mid-section cannot be excited by a 
central load, these higher frequencies are successively 9p,, 25p,, 49p,, 81p,, 
and so on. 


Our purpose now is to follow the steps in the formal argument of Section 2 
et seq. in terms of the concrete example before us; and this will involve certain 
procedures that are necessary only in the standard method of approach, but have 
to be included in order to demonstrate the superiority of the newer one. Referring 
to Fig. 4 (a), we shall suppose the length of the beam BC to be 30 ft., and that the 
central load Q is evenly distributed over the middle foot, which means that the load 
per foot over the loaded length is 


The first step is to express the loading distribution P (x) of Fig. 4 (a) as a series 
of “ normal” loadings L, (x), L, (x), L; (x), ...-., L, (x), ......°0, aS indicated in 
equation (1) but with even values of m omitted. It can be shown that in the 
present example 


so that L, (x), L, (x), ......, L, (x), ...... are given respectively by the loading 


diagrams of Fig. 4(b). As mentioned in Section 3 it next becomes necessary in the 
standard method to decide how far along this series of loadings we have to proceed 
before we are justified in assuming that the sums of the terms give a loading that 
is a near enough approximation to the true loading of Fig. 4(a). To judge this, 
we examine Fig. 4 (c), where the result of adding successive terms up to the fourth 
is shown as loadings (A), (B), (C) and (D) respectively. It is seen that, as an 
adequate representation of the true distribution of 4 (a) for the purpose of obtaining 
shears and bending moments, curve (D) obtained from the first four terms of the 
series is wholly unsatisfactory: a good many more terms have to be included before 
anything like a close approximation is reached. This is the point made in Section 3. 
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Suppose, however, that in order to operate the standard method a few more 
loadings are added to the series shown in 4(+). The next and last step is to 
determine separately the dynamic effect of each of these loadings, L, (x) sin p,t, 
L, (x) sin p.t, before adding them together to obtain 
the resultant. 


Such a roundabout procedure is clearly wasteful when, as often happens, only 
one or two of the lowest modes have any appreciable dynamic effect, 
i.e. inertia forces. 


In the present case the first normal loading L, (x) sinpt of Fig. 4(b) excites 
inertia forces with the same spatial distribution, but varying in magnitude and in 
time according to the shaded area of Fig. S(a). It is seen that at time 7, the 
combined effect of applied and inertia loading is greater by some 25 per cent. than 
the maximum value of the applied first normal loading alone. In the next mode 
Fig. 5 (b) shows that the inertia loading excited by L, (x) sin p.t is negligibly small. 
For all higher modes, therefore, the normal loadings L, (x) sin p.t, L, (x) sin p.t, 
etc., excite no inertia forces, so that the applied loading is resisted in each case 
by the elastic forces of the system alone. 


The total effect of the transient load Qsin p,t may thus be exhibited (cf. 
equations (2) and (3)) as the sum of the following sets of forces 


1. =L, (x)sin p.t+(1.F.), [varying with time as in Fig. 5 (a) ] 
2. =L, (x)sin p.t+(negligible inertia forces) 


3. =L,(x)sin p.t+(negligible inertia forces) 


or (1+2+...... sin p.t+(LF.),. . » 


In other words, the effect is the same as if the load Q were a static load varying with 
time in accordance with the time function sin p.f, except that the inertia forces excited 
in the first mode (and in this case in that mode alone) have to be added as extra 
static forces varying with time, as shown in Fig. 5 (a). 


By following this procedure it is seen that the necessity for considering the 
higher mode normal loadings represented by L, (x), L; (x), etc., is obviated and— 
of comparable importance—that the .error involved in representing the actual 
distribution of the applied load by a limited number of the higher mode loadings 
is also avoided. 
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J 
— 
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Fs} = CRITERION OF STRESS 
a 
fe) 
(I.F.), INERTIA LOADING VARIATION, WITH TIME, 
2 (SHADED) OF APPLIED LOADING 
< 
< 
a. 
= 
LLL’ 
2P. Pp P 
(a). 
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Fig. 5. 


Showing how dynamic (i.e. inertia) effect is important only in first mode and its negligible value 
in the second mode. 


(a) Time variation of applied and of inertia loading in first mode (both spatially distributed as 
Li(x) of Fig. 4(d)). 


(b) Time variation of applied and of inertia loading in second mode (both spatially distributed 
as L(x) of Fig. 4(d)). 
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7. SOME SHORT CUTS IN PRACTICAL PROBLEMS 


In actual practice it may often be convenient to expand the time function F (1) 
into a Fourier series giving 


F (t)= A, sin (int /T) 
i=l 


where T is the total duration of the transient force, i is a general symbol and r is as 
small as possible consistent with the sum of r terms being a fairly good representation 
of the function F (ft). 


The integral in (6) then becomes 


t 


{ cos Pn(t—t,) & (ix/T) A, cos (int, /T) dt, . 
i=l 


which is easy to deal with analytically. 


It is of interest to note a possible short cut in the process of representing F (rf) 
by a Fourier expansion, and the same short cut applies to the use of F (f) in its 
unexpanded form. This consists of dividing the function F (ft) into two or more 
distinct parts as shown in Figs. 6 (a) and 6 (b) where the original curve is divided 
into two parts. 


Figure 6 (a) shows the curve of F (t) as a function of ¢, having a moderately 
even contour, apart from a sharp perturbation over the region CD. This curve can 
be divided into two components, one F, (t), which is everywhere smooth and has 
the same duration as the original F (t), and another F, (t) which has a much smaller 
duration. It requires no mathematical proof to see that the representation of 
curves (b) by two separate Fourier series over different time intervals is a far easier 
task than representing curve (a) by a single series to the same degree of accuracy. 


A parallel computational economy is achieved if F,(t) and F,(t) are used 
in their unexpanded form. For F, (ft) is smooth enough not to excite more than 
a few of the graver modes of the system, while F, (t) is sharp enough to excite a 
few of the higher modes without affecting the graver modes. This will be the case 
particularly if the total impulse of F, (t) is zero, i.e. if the area under the curve 
F,, (t) relative to the base line is zero. 


Note. A most useful short cut is made by realising that it is never necessary 
in practice to compute the component normal loadings L, (x), L, (x), etc., for the 
purpose of estimating their inertia effects: it is only necessary to determine the 
shape of the corresponding normal modes of displacement. Thus it is not necessary 
to calculate the value of L, (x) before we can find the bending moments it produces 
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F(t) 


{a) (b) 
Fig. 6. 


and hence, the bending moments of the corresponding inertia forces (I.F.),. This 
follows from the fact that, so long as the beam (or system) is constrained to deform 
in the n' mode, the bending moments produced by the n™ normal component 
loading of the applied load P (x) are identical with those produced by the applied 
load itself. 


Suppose for example, that the curve in Fig. 7 (which must be obtained either 
by calculation or from a resonance test) represents the shape of the first (or funda- 
mental) normal mode of the beam AB having a mass distribution m (x). We require 
the bending moments produced in this mode by the first normal loading L, (x) 
without computing the value of L, (x) (which resulted in equation (10) in the 
numerical example above). All we have to do is to equate the strain energy in 
the beam to the work done by Q, i.e. Qy,/2, from which y, (the displacement of Q) 
is determined, and hence the value of y, and so the bending moment, at all 
other sections. 


In calculating the internal energy, moreover, it may be noted that it is far 
l 


more accurate and convenient to compute the kinetic energy [ m (x) y*dx rather 


than the strain energy > (33 dx to which it is equal. 


PART II 


8. APPLICATION OF THE METHOD TO IRREGULAR SYSTEMS 


The segregation of the inertia forces from the total resistance to the applied 
transient force, which is the basis of the above method, depends upon finding for 
each mode the Equivalent Static Force (E.S.F.) and expressing the Inertia Force 
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Fig. 7. 


as the difference between the Applied Force and the Equivalent Static Force. The 
applied force is known, being part of the data of the problem, and the E.S.F. for 
each mode is expressible in terms of the frequency alone. However irregular the 
elastic system may be, its frequencies in the various modes, if not calculable, can 
always be found by experiment. It follows, therefore, that the inertia forces can 
be expressed as certain static forces in the various modes. 


The difficulty with the irregular system is in translating the static forces in the 
various modes into displacements and stresses in the system, and the main source 
of that difficulty is that we cannot define the modes analytically, nor measure them 
experimentally. This point will be made clear by considering a practical example 
such as an aeroplane with a heavy gun mounted in the fuselage nose. Certain high 
frequency modes are excited in the structure adjacent to the gun mounting by the 
transient forces from the gun reaction. The frequencies of these modes can be 
found experimentally by applying a harmonically oscillating force along the gun 
barrel at gradually increasing frequencies and noting the resonance frequencies. 
The modes themselves cannot be found nor, therefore, the relation between the 
displacements in any mode and the static force at the gun necessary to produce them. 


The engineer is not so much interested in displacements as in stresses and it is 
possible to relate the stresses in any mode directly to the static load. The 
procedure is as follows :— 


At, and in the neighbourhood of, each resonance frequency, the amplitude at 
the point concerned—in this case at, and in the direction of, the gun barrel—is 
noted; the magnitude of the oscillating force and the stresses at various relevant 
points in the surrounding structure are also noted. The latter would be measured 
in practice by the use of a number of resistance strain gauges in conjunction with 
an oscillograph. 
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The purpose of measuring the amplitude at neighbouring frequencies is to 
obtain the shape of the resonance peak as accurately as possible so as to deduce 
from this the amount of damping present at the particular resonance frequency. 


The differential equation for a damped system under forced vibration may 
be written in the form 
mz+cz+kz=F, sin of 


and the steady amplitude is 


F,/k 


where k= the force per unit static displacement along z. 
w, = resonance frequency for the n‘* mode. 
c. = critical damping. 
By taking two values of z corresponding to two values of w in the same resonance 
peak, k can be eliminated and the value of c/c, obtained. 
We have then, at reson._.ice 
from which k may be found, since c/c. is now known as well as the oscillating 
force F,. 


Thus the relation between stresses and unit static displacement in any mode 
is established, enabling us to find in turn the stresses corresponding to the inertia 
forces in the various modes, and hence to apply the method given in Part I. 


REFERENCE 


1. WiutitaMs, D., and Jones, R. P. N. (1948). Dynamic Loads in Aeroplanes Under Given 
Impulsive Loads with Particular Reference to Landing and Gust Loads on a Large Flying 
Boat. R. & M. 2221, 1948. 


NOTATION 


P (x) spatial distribution of applied load over the structure concerned at 
maximum value (when F (t)= 1). 
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time variation of magnitude of P(x), (whose curve of spatial dis- 
tribution has a shape independent of time). P (x). F (¢) isa maximum 
when F (t)=1. 

n normal loading component of P(x), subject to the same time 
variation F (t) as P (x) itself. 

time variation of inertia loading in n mode (different for each mode). 
circular frequency of structural system in n'* mode. 

circular frequency of harmonical force Q sin p.f in the example. 
current time as distinct from a particular instant defined by t. 
Equivalent Static Force in the n mode, denoting the steady static 
force (having the same distribution as L, (x) ) that produces the same 
displacement as that actually produced dynamically at any given 
instant in the n'" mode. It is equal and opposite to the elastic 
resisting force. 

L(x) F 

inertia loading in the n'" mode having the same spatial distribution as 
L, (x) but its own time variation f, (¢). 

dF (t)/dt. 

spatial distribution P(x) of applied load in numerical example. 
length of beam. 

general symbol denoting the i’ Fourier component defining F (t) 
total duration of transient force. 

mass. 

damping coefficient. 

critical damping coefficient. 

amplitude of harmonic force F, sin wf. 

elastic coefficient. 

circular frequency of forced vibration. 

resonance frequency of system in n“ mode. 

moment of inertia. 

Young’s modulus. 
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SKIN FRICTION IN THE LAMINAR BOUNDARY 
LAYER IN COMPRESSIBLE FLOW 
by 


A. D. YOUNG, M.A., A.F.R.Ae.S. 
(Department of Aerodynamics, College of Aeronautics, Cranfield) 


SUMMARY 

From an analysis of the work of Crocco and others, semi-empirical formule 
are derived for the skin friction on a flat plate at zero incidence with a laminar 
boundary layer. These formule are 


w-1 
R,=0.664 [0.45 + 0.557 (0) + 0.09 (y — 1) M,?o*] 
for the general case of heat transfer, and 
R,=0.664 [1+ 0.365(y— 1) 2 


when there is no heat transfer. 


The problem of heat transfer and the effect of radiation are discussed in the 
light of these formule. The second formula is then utilised in the development of 
an approximate method for solving the momentum equation of the boundary layer 
on a cylinder without heat transfer. The method indicates that with increase of 
Mach number there is a marked forward movement of separation from a flat plate 
in the presence of a constant adverse velocity gradient. 


NOTATION 


The principal symbols used are listed below. Others are defined where they 
occur in the text. 


Suffix 1 refers to quantities measured at outer edge of boundary layer, suffix 
“w” to quantities measured at the wall, suffix 0 to some standard condition, e.g. 


main stream at leading edge of plate, suffix “i ” refers to incompressible flow. 


Paper received November 1948. 
{The Aeronautical Quarterly, Vol. 1, August 1949] 
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distance measured along surface. 
distance measured normal to surface. 
velocity component in x direction. 
velocity component in y direction. 
pressure. 

standard length (e.g. length of plate). 
Mach number. 

temperature (absolute). 
mechanical equivalent of heat. 
enthalpy (=Jc,T, for a perfect gas with constant specific heats). 
specific heat at constant pressure. 
specific heat at constant volume. 

viscous stress. 

density. 

coefficient of viscosity. 

constant x y/ / (2x). 

coefficient of heat conduction. 
Prandtl number = uc, /k. 

defined by » oc T*. 

2re/p,st,*. 

u,xX/ Vv). 

u,L/v,. 

1/1,. 

u/uy,. 

Py. 

function defined in equation (6). 


function defined in equation (48). 
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g(M,,,o0) function defined in equation (51). 


F (n) function defined in equation (9). 
G (n) function defined in equation (17). 
Q rate of transfer of heat from surface of unit breadth and length L 
to gas. 
q local rate of heat transfer per unit area. 
Sy total frictional force on surface of unit breadth and length L. 
Cz Stefan-Boltzmann’s radiation constant x emissivity of surface. 
x 
6 momentum thickness = (1 )dy. 
piu, u, 
x 
displacement thickness (1 = )dy. 
0 
H / 0. 
v 
0 
a" the value of Y defining the outer edge of the boundary layer. 
p* P| 
6/L. 
u* u/u. 
x* w/ 
A : 
fho 
INTRODUCTION 


The general analysis of the flow in the laminar boundary layer at high speeds 
is of formidable complexity. This complexity derives from the variation of viscosity, 
heat conductivity and density with temperature and the consequent inter-dependence 
of the equations of motion and energy. The simplest problem, namely uniform 
flow past a flat plate at zero incidence, has naturally attracted most attention; 
nevertheless results for particular cases have only been evaluated after long and 
laborious calculations. We are indebted for a store of such results to the work of 
Busemann('), Karman and Tsien’), Hantzsche and Wendt”), Emmons and 
Brainerd‘), Cope and Hartree’), and Crocco). The work of the latter is 
particularly important since, in addition to results for particular cases, it reveals 
results of a general character. The more general problem of the boundary layer on 
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a cylinder is still in large measure unsolved. Valuable pioneer work on line 
analogous to Howarth’s series solution for incompressible flow‘) has been begun 
by Cope and Hartree, and Howarth*) has similarly developed solutions for the case 
when o=1.0, and w=1.0. In addition, approximate methods have been developed 
by Frankl), Oswatitsch and Wieghardt’*, Young and Winterbottom'”), 
Illingworth”), and Howarth), in which solutions are obtained of the momentum 
and energy equations on lines analogous to Pohlhausen’s method"*), or in which 
other assumptions or approximations suggested by incompressible flow 
theory are made. 


In the following a general semi-empirical formula is developed for the skin 
friction on a flat plate at zero incidence, both with and without heat transfer in 
a uniform flow of a perfect gas with constant specific heats. This formula is derived 
from an approximate solution of Crocco’s integral equation for the viscous stress. 
The form of this solution is retained but the constants are adjusted to give the best 
overall agreement with the available calculated results for particular cases. The 
problem of heat transfer and the effect of radiation on the temperature measured 
by a thermometer are discussed briefly in the light of this formula. 


The formula is then applied to the development of an approximate method for 
solving the momentum equation of the boundary laver on a cylinder with zero heat 
transfer. The accuracy of the method cannot yet be gauged, as there are no accurate 
results available with which to compare it. It has two features to commend it; 
firstly, it gives the correct result for the case of zero external pressure gradient; 
secondly, it is relatively simple and quick to apply, the main operation in any given 
case being a simple integration. 


The method has been applied to estimate the distance to separate from a flat 
plate in the presence of a constant negative velocity gradient when o=1.0, and 
w=1.0. This is the case considered by Howarth") by a method closely analogous 
to Pohlhausen’s method. It is doubtful whether either method can be relied on to 
give the separation distance accurately, and their results differ considerably for high 
Mach numbers, but both methods agree in predicting a rapid forward movement of 
separation with increase of Mach number at the leading edge for constant ratio of 
(velocity fall per unit length of plate) to (initial velocity). 


2. SKIN FRICTION ON A FLAT PLATE AT ZERO INCIDENCE 
2.1. Boundary layer equations 
The equations of motion for the boundary 20) are 
du 
pu ox +p 


dy 


(1) 


0. 
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The equation of continuity is 


0 0 
ax (pu) + ay (pv)= 0, (2) 
and the equation of energy is 


when /J=Jc,T (enthalpy), and o is the Prandtl number=yc,/k and is 
assumed constant. 


The equation of state for a perfect gas is 
where R is the gas constant and m is the molecular weight of the gas. Since p 
is constant for the flat plate at zero incidence, p=p (J). 


The viscosity coefficient » is a function of T only, the relation being given by 
Sutherland’s formula 


1.5 


, where c is a constant and is equal to 110°K, for air. 


It is convenient to approximate to this relation by writing 


where w varies according to the range of T in which we are interested, but can be 
taken equal to 8/9 for T between 75°K and 300°K"*). 


2.2. Crocco’s Transformation” 


The case of the flat plate at zero incidence is the only one which permits the 
reduction of the boundary layer equations to ordinary equations in terms of a singie 
independent variable. This variable is of the form £=constant x y/./(2x). In this 
case we can write ]=/ (u), p=p (u), n=, (u) and 


du constant (u) 


If we transform from £ to u as independent variable, the equations (1), (2) and (3) 
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reduce to two equations in f and /, namely 


.-. . « 


where dashes denote differentiation with respect to u. 
These equations can be put into a non-dimensional form by writing 


u | v( 2 ) 
— =n, =i, — =r, — =m, F(y)= =)f@, 9) 
uy Pi Pi fu) ( 


1 


when suffix | refers to quantities in the main stream. 
Equations (7) and (8) then become 


and we note that u,*/1,=(y—1)M,*; dashes now indicate differentiation with 
respect to 7. 


The boundary conditions are 
9=0, F’=0, i=1,/1,; 
g=1, F=0, i=1. 
Suffix “ w ” denotes quantities measured at the wall. 


The equation of state (4) gives 


and the viscosity—temperature relation (5) becomes 
Using (12) and (13), equation (10) can be written 
From (6) and (9) it is easy to deduce that 
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where c; is the local skin friction coefficient (=27,/p,u,*) at a point distant x from 
the leading edge, and R,=u,x/¥,. 


Once F (»)) is determined, the relation between y and u is given by (6) and hence 


| 
0 
which reduces to 


2.3. Hantzsche and Wendt's Transformation” 


Hantzsche and Wendt arrived at different forms of these equations by using 
different non-dimensional expressions involving / and f(u). Thus, they wrote 


j=1/u,?, and G V . . (7) 


and their resulting equations which correspond to (14) and (11) were 
GG” + 2n (jw/ =9, (18) 
and G Gj =0. (19) 


In this case the skin friction coefficient is given by 
1 1 
R.=G (0)(j,/jn) 2 =GO)(T,/Te) 2. ‘ (20) 


2.4. Some important results of Crocco’s analysts 


It will immediately be apparent from the form of equations (14) and (11) (or 
from equations (18) and (19)) that considerable simplifications result when 
pu=constant and hence »=1.0, or when o=1.0. 


Thus, when o = 1.0, equation (11) reduces to 


and hence i is a quadratic function of », independent of F; taking account of the 
boundary conditions this function may be written 
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(22) 
ot i=1(0)+ [1-1 @+ 
Having i as a function of u, equation (14) for F takes the form 

FF” +h(n)=0 


where h(n) is known. For details of an iterative process that can be used to solve 
this equation for F see Ref. 6. 


When pu=constant or »=1.0, equation (14) simplifies to 
which can be transformed to the Blasius equation for incompressible flow, namely, 


ay 


dg +S qe 


q 


by writing F=}4 n= 


The relation between € and y,, the lateral ordinate in incompressible flow, is 


/(@)(3,)- 


From equation (16), the lateral ordinate in compressible flow corresponding to a 
given value of 7 is 


= = nd. = id i 
= | dy, = | a ‘ . (24) 


Hence for this case all solutions are simple transformations of the incompressible 
flow solution, the relation between corresponding values of y for given values of u 
(or n) being given by equation (24). 


We note immediately that since »=0 when £=0, F (0) is the same for 
compressible as for incompressible flow and hence in all cases when pu =constant 
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In the general case, when neither » nor o is necessarily unity, let us suppose we 
have solved for F as a function of 7. Then equation (11) which is linear in 7 can 
be integrated twice to give 


i (n)=i(0)+7 (0) A (y, 7) (u,?/T,) ©), 
where A (n, dn,, 
0 (26) 
" 
and B(y, 7)= Fay F (0) dn, ¢ dn,. 
0 0 
Alternatively, taking account of the fact that /(1)=1, we have 
i(n)=1-7 O)C(y, 7) +e (u,?/1,) DQ, ©), 
1 
where C (n, 7) = dn,, 
(27) 
1 " 
and D(n, Ee dn, dn,. 
q 0 
Hence, for the case of zero heat transfer when 7’ (0)=0 
I@-!10) 
(28) 
and in particular =2cD (0, c). 
u,?/2 


Since D (0, 7) does not involve M,, it will be independent of Mach number for 
all cases where F(n) is independent of Mach number, i.e. for all cases where 
pu=constant (or w=1.0). Pohlhausen’s solution”) for incompressible flow gave to 
a very close order of accuracy 

1()-7(1) | 


u,?/2 
hence 2cD (0, 
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and therefore for all cases where pu =constant, 


However, by a laborious trial and error method Crocco“? integrated equations 
(11) and (14) step by step for the cases e=0.725, a range of Mach numbers up to 
5.0 and values of » of 1.25, 0.75 and 0.5. In every case the results showed the 
functions A, B, C and D to be practically independent of the particular value of « 
assumed, and hence independent of the function F ()) within the range of values 
of » considered. Crocco suggested that the reason for this lay in the fact that while 
the variation of F (y) with » may be considerable, it occurs in the integrals from 
which A, B, C and D are derived in the form Fan or [F oy. The ratio 


F (0) F (0) 
F (n)/F (0) must in every case be unity when »=0, and be zero when »=1.0, and 


hence its possible variation with is limited. Further, as long as o is not far off 


unity, the exponents (1 —) or (o — 1) are small and reduce still further the sensitivity 
of the expressions A, B, C and D to variations in ». The final result is therefore 
that in general the variation of enthalpy with velocity is practically independent of 
the viscosity—temperature relation for values of o not very different from unity, and 
is given by equations (26) or (27) with the functions A, B, C and D calculated for 
the case »=1.0, when F (yn) is given immediately by the Blasius incompressible 
flow solution. In particular, when there is zero heat transfer, equation (29) for the 
wall temperature is of general validity, a fact that is confirmed by the cases calculated 


by Brainerd and Emmons"). 


Since i (n) can now be regarded as independent of F (n) in general, and is readily 
determined from the case »=1.0, the solution of the main equations (10) and (11) 


reduces to the solution of the equation 


FF” +h (n)=0, 


where h(n) is known. This enabled Crocco to develop a more rapid and more 
accurate process of solution than the first process he had used. For further details 
see Ref. 6. 
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2.5. Skin fnction dtistnbution. Deduction from Hantzsche and 
Wendt's results for case of zero heat transfer 
Hantzsche and Wendt’s calculations’? for the case of zero heat transfer showed 
the interesting result that the quantity G (0) was only a slightly varying function of 
» and M, for ~=1.0, and was independent of + and M, when »=1.0. In general, 
therefore, we may expect G (0) to vary little with o, » and M,. The importance of 
this may be seen from equation (20) which with equation (29) gives us 


A reasonable approximation to G (0), reproducing its small variation with », 
and M, should then give us a reliable formula for c,/ R., for which the major 
variation with w, o and M, is described by the denominator of (30). In Fig. 1, G (0) 
is shown as a function of M, for various values of », for the case «= 1.0, and this 
can be represented fairly closely by 


1+ 


G (0) = 0.664 | 


Chis suggests that in general 
5 I = 
1) M,?c* 
and hence for the case of zero heat transfer we may expect from (30) that 


F (0) 


2.6. Skin friction distribution with or without heat transfer. 
Deduction from approximate solution of integral equation 
for F 


From equation (14) and the boundary conditions for F it follows that 


" 


From a survey of Crocco’s results for F it appears that in all cases a very close 
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\ 
ASYMPTOTIC VALUE d 
O's FOR M,=0o 
a 
i-w 
(¥-1) 
FOR o =!1:0 
= 7 
0-2 NO HEAT TRANSFER 
O:! 
2-0 4:0 6-0 8-0 10-0 
M, 
Fig. 1. 
Skin friction of flat plate at zero incidence (no heat transfer). 
The function G (0) according to calculations of Hantzsche and Wendt (c=1.0). 
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fit to F is given by 


(0) 
Hence, from (32) we may expect to get an even closer fit from the relation 


1 


" 
In particular, 
1 " 


We have ij as a function of » and o given by (26) or (27) so that in theory we can 
determine F (0) from (33) in any given case by integration (which must in general 
be numerical), and hence from (15) we can obtain the skin friction coefficient. We 
are interested, however, in obtaining a general explicit formula for F (0) in terms 
of i(0) (or T,/T,), M,, « and ». Analytic integration of (33) appears to be 
impossible in general; the process adopted, therefore, has been to develop a crude 
approximate solution of (33) to suggest the form of the required formula for F (0); 
the constants in this formula have then been adjusted to give the best fit with the 
numerical results of Crocco and other workers. 


Consider the case when ~=1.0. Then we can write 
=A,+A.+A3n’, 
where A,=i (0), A,=1—i(0)+u,?/21,, —u,?/21,. 
Write »=sin x, then (33) can be written 
OM 
F (Oy = [| 2 sin x (A, +A, sin x +A, sin? dy | dy,. 
0 


0 


Let A,=1+i(0)+u,7/2/,, then A, > A, +A, sinx+A, sin’ y. 


Hence, we can write 


2 
0 


[fo sin [1+ (42 sin sin? J] dx,, 


0 


+ 34. “6-3 +terms in (34) 


F 
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This leads to the result that when w= 1.0, 


F (0 =2- = 


N 


and therefore F (0)=0.655 (compare the accurate value of 0.664) . 


Write equation (34) in the form 

F@ _,: [5 A, _ A, 

A, 
A, (32-97) 


124, (4-2) ] + terms in (w — 


Then assuming » is not very far from unity (as will be the case in practice), the 
terms in (w—1) can be allowed for to some extent by approximating to the right 
hand side with 


F (0) 3n—8 As 


3x 20 — 62 u,? [ 217-64 7) 


[0.553 + 0.447 (0)+ 0.096 u,7/1,J2 . 


Even before any attempt is made to modify the constants in this formula to improve 
agreement with calculated results we may note that it gives good agreement with 
formula (31) for the case of zero heat transfer and «7=1.0. Thus, in this case, 
from (29), 


i(0)=1+ =14u,?/21,, 


and hence (35) becomes 


FO). =(1+0.32 u,? 2 [1 +0.32 1) 15 


This comparison indicates that any modification to (35) should retain the 
following features : — 


(i) The sum of the constant term and the coefficient of i (0) should be unity. 


(ii) Half the coefficient of i(0) plus the coefficient of u,?/7, should be 
approximately equal to 0.3. 


The comparison also indicates that to take account of values of « other than unity 
the coefficient of u,*/J, should include the term c*. 
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Formula (35) with this last modification was in fact found to give results in 
surprisingly close agreement with the calculations of Refs. 2 to 6, bearing in mind 
the approximate nature of the analysis. The closest overall agreement was found, 
however, with the formula 


FQ) 
FO... [0.45 +0.55 (0)+0.09 (y—1)M,?o*] 2 . 
Using the correct basic value of F (0)..., = 0.664, this can be written 
R,=0.664 [0.45 + 0.557 (0) + 0.09 (y 1) M205] 2 ; . (36) 


For the case of zero heat transfer this reduces to 
R,=0.664 [1 +0.365 (y—1)M,%o"] 2. . 


Table I compares the results given by equation (36) and Crocco’s calculations, and 
it will be seen that in every case the agreement is within one per cent., which is 
within the accuracy of Crocco’s calculations. In Table II Crocco’s results for zero 
heat transfer are compared with the results given by equation (37), and the agreement 
will be seen to be equally good. Finally, in Fig. 2 some of the results of Refs. 2, 4 
and 5 are compared with equation (37) for the case of zero heat transfer, and in 
every case close agreement is found. One can conclude that the formule (36) and 
(37) are reliable to within about one per cent. for values of the Mach number up 
to 10, and for values of » and o likely to be of practical interest when air is the 
working fluid.t 


3. HEAT TRANSFER AND RADIATION—FLAT PLATE AT ZERO 
INCIDENCE 


3.1. Heat transfer 

If Q is the total heat transferred per unit time from a surface of unit breadth 
and of length L, say, g is the local rate of heat transfer per unit area and Sy is the 
total frictional force on the surface, then it readily follows from Crocco’s analysis 


that for o near unity 


Q q_ se 
Se 2/8 


where 7,, is the temperature measured by a thermometer, i.e. 7,,=—7, for zero 


+ The coefficient of (y—1) Mc} in equation (37) is somewhat greater than that suggested by 
equation (31) based on Hantzsche and Wendt’s results. However, the complexity of the 
computations involved leaves it likely that small systematic errors may have been present 
in the latter's calculations. 
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Comparison of Crocco’s calculated values of skin friction on a flat plate with values 
given by equation (36) (¢ =0.725). 


TABLE I 


| (Crocco) (Equation (6) 
pas tel 125 | 025 0.940 0.942 
| 2 125 «| 025 0.957 0.960 
| 5 1.25 | 0.25 1.040 | 1.039 
| 1 | 1.006 1,004 
2 125 10 1.016 LOIS 
| 5 125 | 10 1.076 | 1.074 
| 1 125 | 20 1.056 | 1.059 
| 2 125 | 20 1.066 | 1.066 
| 5 125 | 20 1.11 | L111 
pes 0.75 | 0.25 1.070 | 1.062 
| 2 0.75 | 0.25 1.049 1041 
| 5 0.75 | 0.25 0.960 | 0,963 
oie 0.75 1.0 0.985 | 0,986 
5 0.928 0.931 
1 0.75 | 20 0.946 0.944 
2 0.75 2.0 0.940 0.938 
5 0.75 2.0 0.903 0.900 
1 0.5 0.25 1.139 1.128 
2 0.5 0.25 1.098 1.089 
5 0.5 0.25 0.931 0.927 
1 0.5 1.0 0.991 0.992 
2 0.5 | 1.0 0.970 0.971 
5 oe +e 0.868 0.868 
| 1 0.5 | 2.0 0.897 0.892 
| 2 0.5 |} 20 0.886 0.879 
| 5 0.5 | 2.0 0.815 0.811 


+ To minimise the effects of a possible error of the order of one per cent. in Crocco’s 
calculations, ratios rather than absolute values of skin friction coefficients are compared. 
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heat transfer =T, E + | 


If we now use equation (36) to derive Sp and 7,, we find 


Sy = 0.664 [0.45 + 0.55 i (0)+ 0.09 (y—1) M,?o*] . G9) 

and T~=0.332 [0.45 +0.55 i (0)+0.09 (y—1)M,?o%] 2. . (40) 


Equations (38), (39) and (40) enable us to calculate g and Q in any case where 
,, M,, and » are given. 


3.2. Effect of radiation on thermometer temperature 


The following treatment is a modification of that due to Hantzsche and Wendt”), 
use being made of the explicit formula for F (0), equation (36). 


Allowing for the radiation of heat in the thermometer problem, then the rate of 


TABLE Il 


Comparison of Crocco’s calculated values of skin friction on a flat plate with values 
given by equation (37) for case of zero heat transfer (c =0.725). 


M, | | OW Rs _ | Re 
| | (Equation (37) 
1 1.015 
2 | 125 | ~ 1053 1.052 
5 1.193 
1 | O75 | 0984 0.986 
2 | 075 | 0950 0.951 
5 | o7 | 0942 0.838 
1 os 0.969 0.971 
2 | os | 0908 0.904 
| 5 05 | 0.707 0.702 
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heat transfer from gas to wall must equal the rate at which heat is radiated. Hence 
L 


Q= Cz (T,‘ T,*}ax, 


0 
where Cy, =Stefan-Boltzmann’s constant multiplied by the emissivity. 


Using equation (38) we have therefore 


P Tw T, n 2F 0 


(0)—1— M204] 


or 


Then K,= —— ~ F (0), 


(42) 
where, as before, we have 


F (0)=0.664 [0.45 +0.55 i (0) +0.09M,204]- 


The problem is, in general, given Kz, what is the value of i(0)? With the aid of 
equation (46) the problem is best treated in reverse; that is, for given values of M,, 
w, « and i (0) the corresponding values of Kz can be calculated and the values of i (0) 
can be plotted against Ky. We may note immediately that for Kz=0 (as when 
L=0), Tx =T, for Ky (as when L=00), i(0)=1, i.e. T,=T,. For intermediate 
values, i (0) decreases steadily with increase of Kp. 


4. THE LAMINAR BOUNDARY LAYER ON A CYLINDER 


The method to be developed here is essentially the same as a method previously 
put forward by the author"), but with a number of improvements suggested by the 
results given in Section 2. 


Suffix 0 will be used to denote undisturbed stream values, and suffix 1 to 
denote values just outside the boundary layer. 


The momentum equation of the boundary layer is readily shown to be"? 
u, Pi 
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where (1 dy, 8 ) ay, 


H =8**/6, and dashes denote differentiation with respect to x, the distance along 
the surface. 

We will now follow some of the implications of an approach on the lines of 
Pohlhausen’s method for incompressible flow. Instead of expressing the velocity 


distribution as a function of y, however, it will be expressed as a function of 


Y= { - ~ since we know that on a flat plate, with » = 1.0, the velocity distribution 


expressed as a function of Y is independent of Mach number. Thus, suppose the 


velocity distribution across the boundary layer to be given by the quartic 


with the boundary conditions that 


for Y=3,°, = = =a =(); 


oT Ou 
and for Y=0, u=0, oY =0, and hence 


The quantity 5,* is at present unknown, but is assumed to correspond to the value 
of Y defining the outer edge of the boundary layer. 


The first boundary layer equation yields 


But ay u ay = [lw since (34 
O7u 
Uw 


With these boundary conditions we can solve for the coefficients a, b, c and d in 
terms of 5,*, and we obtain 
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_4(I2+A) -wA _-u,(4—-A) ,_u,(6—-A) 
(45) 
where 
w 2 A 
tw _ Me (12+A) (46) 


piu,” 6p 


From this point the method departs from the classical Pohlhausen approach. It 
depends on the solution of the momentum equation, making the following 
assumptions 


(i) Equation (46) is accepted. 


(ii) The variation of H with local velocity, pressure gradient, » and o is 
neglected, H being treated as a function of M, only. The justification for 
this lies in the fact that H enters the solution of the momentum equation 
in a form suggesting that the solution is relatively insensitive to small 
variations of H, and for slim bodies at small incidences the possible 
variations of H are unlikely to be large. 


For a flat plate at zero incidence, with » and o both unity we know that u/u, 
is a unique function of Y, and also the total energy is constant across the boundary’ 
layer. Hence, in this case 


where suffix “i” denotes the value in incompressible flow. 


Also =| (1 dy = 

But = T, =1+ M, 1 


+ It is possible to avoid making this assumption, but the computing involved then becomes 
considerably heavier. 
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+ M.? (&:** +4). 
Therefore 


The Blasius solution gives 
H,=2.59, 


This relation is in almost exact agreement with the results of the calculations of 
Brainerd and Emmons“). 


(iii) The ratio 5,*/@ is also assumed to be independent of variations in local 
velocity and pressure gradient. For the present we write 


5,*/0=f (M,, o)t : . (48) 


making no further assumptions at this stage about the form of the function, 
except to note that we expect it to be constant when » = 1.0. 


From (46) and the expression for A in (45) we have 


Hence, the momentum equation (42) becomes 


u,’ p,’ u,’ 
( Pi 6U, Mo f 


t This function should not be confused with the function f of equation (6). 
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Multiply both sides by »,76, and we get 


2p2 4 _f PoP 1 


At this stage a further assumption is made, viz. the value of u»/, in the second 
bracket is taken to be that for a flat plate at zero incidence. This gives 


ays 
= [ 1+ om. | — 
Write H+2- 9 (M,,w,o)/2,say. . . (1) 


Then (49) can be integrated with respect to x to yield 


=. z. 


1 


0 


1 


The leading edge is taken at x=0, where either u, =0, or u, is finite. If u, is finite, 
then 6=0 there, since we cannot have a finite momentum loss there. Hence 


2 4u. g-1 


We may note here that g occurs as an exponent of u, inside the integral, and as an 
exponent of u, (x,) outside the integral in the denominator. This suggests that 
provided the variation of u, with x is small, small errors in g will have little effect, 
and this provides the justification for the neglect of the variation of H and u./u, 
with pressure gradient and local velocity. 


We have yet to determine the function f. 
For the flat plate at zero incidence (53) reduces to 


4uXx, 
uot 


6/R._ | 4 
d 
ox JR. 
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96 
But C= Polte? 2 ox 
Hence Re=0.664 [1 +.0.365 (y - 1) from equation (37). 
It follows that 
f=9.072 [1 + 0.365 (y — 1) ‘ (54) 
The function g is given by 
g=2(H+2)— 
3 Mo 
af f 


It is sometimes convenient to express these relations in non-dimensional forms. 


where L is some standard length, then (53) becomes 
4 
1 1 
0 
where R= p,L/ vo. 
The relation between p,* and u,* is 
1 


Having u,* as a function of x*, we can then obtain 6* at any point from equations 
(54), (55), (56) and (57), the solution of (56) being obtained by graphical or 
numerical integration. 
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To obtain the local skin friction coefficient when 6* is known, we have 


 (A+12)u, 


and A = . (58) 
=R 


For w* = w/o, We May use equation (50) but it is probably more accurate at this 
point to assume that the local relation between 7,, T, and u,, is the same as for 
the flat plate which leads to the relation® 


2 


The accuracy of the method outlined above cannot as yet be assessed; this must 
await the development of a method that can be accepted as accurate and its 
application to a few test cases. In favour of the above method it should be noted 
that it gives the correct answer for the flat plate at zero incidence and the numerical 
work involved in its application is relatively small. It may also be noted that when 
applied to incompressible flow’) the method gave results which were in very close 
agreement with results given by the standard Pohlhausen method. 


5. SEPARATION OF THE BOUNDARY LAYER 


It is unlikely that the above method will be any more successful in predicting 
the separation point of the laminar boundary layer than is the Pohlhausen method 
in incompressible flow. Nevertheless, it should provide a guide as to the effect of 
increase of Mach number on tendency to separation. 


From (45) separation occurs when A = — 12, i.e. when 


Ré**p, = from (58). 


Hence, from (56), separation occurs at x,, if 


_ 
= | p.tut tat. . - (60) 
0 


u,* 


Howarth has considered the case of a uniform adverse velocity gradient with »= 1.0 
and «=1.0, and has estimated the separation distance for different main stream 
Mach numbers at the beginning of the plate, using a method that is essentially an 
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analogue of Pohlhausen’s method. It is of interest to compare his results with the 
results given by the present method. 
With w= 1.0, and c=1.0, 
f=9.072, g=6.156+0.831M.’, and py* = 1+4(y— 1) M.?. 


Suppose u,=u,— 8x, where 6 is a constant. 


Then u,*=1-— x*, and u,*’= — 


=a, Say. 


Equation (60) becomes 


z,° 


9.072 (1+ (1 ax,*} p,* (x) 


and, from (57) 
“a =i 
p,*= [1 + (22x* — 
Let xx* =s, then we have 
(y- 1) 
[ 1+ — M, 
0.331= — X 
(1—s,¥ [1 + Du. (Q2s,-s,2) 
1 
[1+ Qs—s*) (l—s¥-'ds, . . (61) 


0 


with g=6.156 + 0.831M.,’. 


The process is then, given M,, to evaluate the right hand side numerically for a range 
of values of s,, and to interpolate to obtain the value of s, for which (61) is satisfied. 
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The results for the separation distance x, given by this method are compared 
with those given by Howarth’s method in the following table : — 


— Uo 
M. Howarth Present method 
0 0.156 0.165 
] 0.148 0.148 
3.16 0.107 0.081 
10 0.052 0.013 


The agreement between the two methods is poor for the higher Mach numbers, but 
neither method can be claimed a priori as very reliable for predicting the point of 
separation, and it is impossible at this stage to estimate their relative degrees of 
reliability for this purposet. It can be inferred from both sets of results, however, 


that the separation point moves forward with increase of initial Mach number for 
du 

constant 
Uo dx 


plate) to (initial velocity). 


U,, i.e. for constant ratio of (velocity fall per unit length of 
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THE MEASUREMENT OF GAS TURBINE 
COMBUSTION EFFICIENCY BY GAS ANALYSIS 


by 


L. J. RICHARDS AND J. C. STREET 
(Thornton Research Centre, Shell Refining and Marketing Company Limited) 


SUMMARY 

Following a brief discussion of various methods of measuring combustion 
efficiency in gas turbine chambers, the superiority in accuracy of applied gas 
analysis is stated. Factors leading to inefficiencies in combustion are outlined and 
the exhaust constituents requiring measurement are deduced. The chief difficulty 
in the past in the application of gas analysis methods to combustion efficiency 
determination has been the time consumed in actual analysis of samples; at 
Thornton Research Centre means of overcoming this difficulty have been investigated 
and several rapid and accurate analytical techniques have been developed and 
used successfully. These techniques are described, together with other promising 
methods at present under development. The broad requirements of obtaining 
exhaust gas samples for analysis are discussed and typical sampling systems 
described; precautions necessary in applying these systems are given. Formule for 
the calculation of combustion efficiency from analytical results, together with any 
assumptions made, are given in the text and their derivation explained in an 
Appendix. The techniques described have been developed primarily for use in 
conjunction with rig testing of aircraft gas turbine combustion chambers; the broader 
application—for instance, to engine testing—of the principles involved are 
considered briefly. 


1. INTRODUCTION 

It is well known that in the development and testing of gas turbine combustion 
chambers, one of the most important performance criteria requiring evaluation is 
“combustion efficiency”—a measure of the proportion of the total fuel supplied 
usefully employed in giving energy to the working gases of the turbine. 


The variety of methods by which measurement of combustion efficiency has 
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been attempted has resulted in several different conceptions and definitions of this 
term, but since the primary considerations are those of heat or energy, perhaps th« 
one most appropriate, and therefore adopted in this paper, is the ratio of the heat 
actually developed in the chamber to the heat which would be developed by complete 
combustion of all the fuel supplied to the chamber. In the experimental deter- 
mination of combustion efficiency, calculation of the maximum heat or energy 
available is based on the calorific value and rate of supply of fuel to the combustion 
chamber on test, but measurement of the actual heat generated has been attempted 
in several ways, the most important of which are: — 


(@) by combustion chamber thrust measurement 
(ii) by measuring gas temperature rise in the combustion chamber 


(iii) by chemical analysis of the combustion chamber exhaust gases and 
computation of efficiency from the measured quantities of various 
constituents. 


In the combustion chamber thrust measurement, efficiency is determined by 
comparing the actual thrust, obtained under given chamber operating conditions, 
with the theoretical maximum thrust calculated from fuel and air supply rates, 
calorific value of the fuel, and the exhaust nozzle dimensions. Direct comparison 
on these lines will not give values of combustion efficiency in accordance with the 
earlier definition, but will include the aerodynamic efficiency of the chamber. 
Further tests to determine the pressure loss characteristics of the combustion system 
will enable combustion efficiency to be partially separated from the combined term, 
but the validity of utilising “cold” pressure loss data in the correction of experimental 
results obtained with high temperature gases is doubtful. The method then, although 
perhaps of use for tests of comparison, cannot yield results of any great absolute 
accuracy. Probably the most important disadvantages in the application of the 
principle are the mechanical difficulties involved; obviously, for any acceptable 
degree of accuracy small changes in thrust must be readily detectable, and the 
physical proportions of the necessary flexible air ducts, etc., leading to the chamber 
militate against this. It can be said, then, that thrust measurement shows little 
promise as an accurate means of combustion efficiency measurement, and is at 
present rarely used for that purpose. 


By far the most widely used method for the evaluation of combustion efficiency 
is that of measuring the gas temperature rise in the combustion chamber and 
comparing the results with the theoretical value based on complete combustion. 
A direct comparison of such temperature increments results in a value of efficiency 
slightly different from that given by the “energy” definition, because of variations 
with temperature in the specific heats of the gases under study, but in most practical 
cases numerical discrepancies are negligible, particularly at high values of efficiency. 
The practical application of this principle requires the measurement of gas 
temperature at the outlet of the combustion chamber where, unfortunately, the 
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temperature and gas velocity distributions are seldom uniform. Common practice, 
then, is to traverse the outlet of the combustion chamber with a suitable thermo- 
couple and pitot-static instrument and thus obtain a mean value of the outlet gas 
temperature weighted according to gas velocity distribution. The necessity of such 
a procedure is the prime factor in limiting the accuracy of this method; the accurate 
measurement of high gas temperatures in rapidly moving streams is a problem in 
itself, while, because of the steep temperature gradients frequently occurring across 
the chamber outlet the estimated mean value of indicated gas temperature is often 
open to doubt. Further errors must be expected in estimating the maximum gas 
temperature rise available from the fuel supplied, for this must be calculated from 
measured air and fuel supply rates, both of which are limited in accuracy, particularly 
the former. 


So far only the unfavourable characteristics of the “ heat balance ” technique 
have been discussed, but its wide application testifies to its general value. In the 
hands of experienced personnel, the method can be made relatively speedy and 
reliable within known limits of accuracy and, although perhaps not suited for highly 
accurate determination of, or the detection of small changes in, combustion 
efficiency, it is at present the most convenient tool for general combustion 
chamber development. 


Where a high degree of accuracy is required in efficiency determination, exhaust 
gas analysis methods are unsurpassed, being both sensitive and accurate on an 
absolute standard. The principles involved consist broadly in the estimation of 
both the heat actually released and the maximum heat available in the combustion 
chamber from the composition of the exhaust gases, and in some circumstances no 
other measurements are necessary. One of the chief objections in the application 
of these principles, however, has been the time consumed in the analysis of samples, 
and much effort has been expended in recent years in devising rapid analytical 
techniques of the required precision. In this paper some of the more successful 
results of such work at Thornton Research Centre are summarised, and their 
application in the determination of combustion efficiency in aircraft gas turbine 
combustion chambers is described; analytical methods envisaged for future use or 
at present under development are also included. 


2. THE PROBLEM 

In contemporary combustion chambers, because of practical limitations, the 
combustion process is never theoretically complete. This means that a proportion 
of the total available heat is not released; this unavailable heat is known as the 
combustion loss. 

Several factors contribute to the combustion loss, the more important of 
which are: — 

(i) The chilling of the combustion reactions occurring in the primary zone, 

mainly by the admission of the secondary and tertiary air. 
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(ii) Slow or delayed combustion, which is usually caused by bad mixing in the 
primary zone and is often indicated by a long flame. This will cause the 
effect of the chilling to be much more severe, as the reaction will be chilled 
at an earlier stage. 


The above two factors will cause the presence in the exhaust gas o! 
partially burnt products. 


(iii) Failure to ignite some of the fuel or products of the secondary reactions, 
ie. cracking products. This again may be due to imperfect mixing 
conditions in the primary zone. This factor will result in unburnt fuel, as 
well as partially burnt products, appearing in the exhaust gas. 


Taking into account these various mechanisms, the likely constituents of the 
exhaust which contribute to combustion loss are: — 


(a) Hydrocarbons (unburnt fuel, cracked fuel, and so on) 

(b) Carbon monoxide 

(c) Hydrogen 

(d) Other intermediate combustion products (aldehydes, acids, peroxides) 
(e) Carbon. 


The relative importance of these may vary considerably with the combustion 
chamber operating conditions, but generally speaking the first three are by far the 
most important. The amount of carbon formed even under heavy coking conditions 
corresponds to a thermal loss of 0.1 per cent.) or less of the total available heat and 
may be neglected. Peroxides have not been detected in normal combustion. Formal- 
dehyde and formic acid‘*) can be detected by direct chemical methods, but again 
their contribution is normally not significant (< 0.1 per cent.). Even when the 
total loss is large and their contribution rises to about 0.3 per cent., this still forms 
a very small fraction of the total loss and may be neglected. For purposes of 
rapidity and simplicity, therefore, constituents (d) and (e) may be neglected for all 
practical determinations of combustion efficiency. 


In all determinations of combustion efficiency on gas turbine combustion 
chambers, there is a conflict between two factors, i.e. rapidity and accuracy. A 
complete and comprehensive analysis could be made of every possible exhaust gas 
constituent, which would enable the combustion efficiency to be determined with 
very high accuracy, but the time involved would be far too long and this would not 
be practicable as a routine test. Also, the accuracy obtained would probably be 
meaningless because the overall accuracy is controlled by the difficulty of obtaining 
a truly representative sample of the exhaust gases. 


Accordingly, bearing in mind the necessity of maintaining an accuracy commen- 
surate with the requirements of the test, various approximations may be made which 
will enable the determination of efficiency to be made with greater rapidity than 
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otherwise. Thus, as stated above, the loss due to carbon, aldehydes, acids and 
peroxides may be neglected. This means that the exhaust gas must be accurately 
analysed for hydrogen, carbon monoxide and hydrocarbons. The method of 
calculation of combustion efficiency from the analytical results, as well as the various 
approximations which may be used, are shown in Appendix I and discussed in a 
later Section. If the requirements of the test permit, an approximation to efficiency 
may be made by neglecting the hydrogen, which normally is only present in small 
amounts, and by measuring the loss caused by carbon monoxide and hydro- 
carbons alone. 


Finally, if rapidity is of the utmost importance, as is often the case in routine 
work, instead of determining the actual thermal efficiency an estimate of this 
efficiency may be made by determining the CO, combustion efficiency. This may 
be defined as: — 

Percentage by volume of CO, in exhaust gas 


Theoretical percentage by volume of co, for complete combustion to Co, and H,O 


This “efficiency” will normally be lower than the true efficiency because it 
involves treating the intermediate combustion products as a total loss, whereas in 
fact they are not so; for instance, the formation of carbon monoxide is accompanied 
by some heat release. The fact that the presence of any free hydrogen in the exhaust 
gases has been neglected would tend to minimise this discrepancy. The relative 
efficiencies obtained when comparing various fuels and conditions are not greatly 
affected by the above approximations, and the rapidity with which results may be 
obtained is greatly increased since, in the case of the CO, combustion efficiency, 
as will be shown later, all that is required is two determinations on the exhaust gas 
for carbon dioxide. 


3. ANALYTICAL METHODS FOR MEASURING IMPORTANT 
CONSTITUENTS OF LOSS 


This Section represents a brief survey of some possible methods for the 
determination of constituents of the exhaust gas which are important in combustion 
efficiency determination. For consideration of space the methods are not described 
in detail, emphasis being placed mainly on the principles involved. 


3.1. Cambridge vacuum method 


This method is primarily used for the estimation of gases which contain carbon 
in the molecule, although the analysis does include an estimation of oxygen and 
nitrogen.. The layout of the apparatus is shown diagrammatically in Fig. 1. The 
principle of the method is simply, by fractional combustion, to convert in turn 
each carbon-containing constituent or group of constituents to carbon dioxide (CO,), 
to condense the CO, so formed and, after distillation, to measure the amount of 
this CO, by the pressure exerted in a known volume. 
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Fig. 1. 
Vacuum gas analysis apparatus. 


The operations involved in converting each constituent in turn to CO, can be 
summarised briefly as follows : — 
(i) The gas sample is admitted to the evacuated apparatus through a trap 
at — 80°C. 
(ii) Any hydrocarbons condensed in the trap are estimated by distilling, 
circulating over CuO at 800°C. and measuring the CO,; the water formed 
in this and subsequent operations is condensed in a trap at —80°C. 


(iii) The gas pressure is measured in a known volume. 


(iv) The CO, content of the exhaust gas sample is estimated by circulating 
the sample through a liquid-air trap and subsequently measuring the CO, 
after removal and storage of the uncondensed gases. 


(v) Any hydrocarbons which remain uncondensed at —80°C. but which 
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condense with the CO, are estimated by passing the CO, over CuO at 
800°C. and measuring the increase in pressure. 

(vi) The CO in the uncondensed gases from (iv) is then oxidised to CO, over 
CuO at 320-330°C. The hydrogen is oxidised at this stage but is not 
measured. 


(vii) The CH, is similarly estimated by oxidation at 1,000°C. The oxygen 
content can also be determined at this stage by absorption by free copper. 


(viii) The residual gas is nitrogen. 


Advantages 


(a) The fact that the CO, produced by each constituent is measured directly 
by the pressure it exerts in a fixed volume, instead of by the difference 
between two large quantities after absorption, means that small amounts 
can be measured very accurately. 


(b) The absence of liquid absorbents. 


(c) The accuracy of the method can be adjusted to suit the quantity to be 
determined, by alteration of the volume in which the final pressure 
is measured. 


Disadvantages 


(a) Time required. For the complete analysis as described, and using a mercury 
manometer to measure pressure, the time required is 3 to 34 hours. This 
time can be shortened if it is only necessary to determine CO, combustion 
efficiency; the analysis then takes less than 2 hours. 


(b) Omission of hydrogen determination. Hydrogen can be estimated by 
measuring the contraction of a sample after oxidation on a platinum spiral, 
but to measure such a small quantity accurately would require a more 
sensitive pressure measuring device and this would increase the time 
required for a complete analysis to impracticable dimensions. Accordingly, 
in the present work, the hydrogen determination is omitted. 

As normally hydrogen is less than 0.1 per cent. by volume, this 
omission is not serious, but it must be emphasised that it is only permissible 
for certain work. General comparative work will be little affected. As 
will be mentioned later, an alternative method for hydrogen determination 
is being developed which will repair the omission in the vacuum method. 


Accuracy 


The normal accuracy is +0.005 per cent. for any carbon-containing gas, 
+0.01 per cent. for oxygen, and +0.05 per cent. for hydrogen. This corresponds 
to the estimation of CO, combustion efficiencies to +0.2 per cent. at an air/fuel 
ratio of 60/1. 
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Fig. 2(d). 
Modified Cambridge vacuum gas analysis apparatus. Apparatus for measurement of CO. 
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3.2. Modified Cambridge vacuum methods 

As mentioned above, one of the chief disadvantages of the Cambridge vacuum 
method is the time required even for the determination of CO, combustion efficiency. 
Accordingly, special attention has been paid to modifying the original method so 
as to reduce the time while retaining the accuracy. 


The method finally adopted) is to draw the gas by means of a vacuum pump 
from a sample bulb of known volume through a trap cooled in liquid air (see Fig. 2). 
The rate of evacuation is carefully controlled by means of a capillary restriction. 
After complete evacuation of the system the trap containing the CO, is transferred 
to another apparatus where the CO, is vaporised and the pressure it exerts in a 
fixed volume is measured. By this means the time-consuming circulation processes 
of the original vacuum method are eliminated. 
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The time for a single analysis is about 30 minutes, but when a series of analyses 
is performed together the average time for one analysis is about 20 minutes. The 
accuracy of the method is at least as good as that of the Cambridge vacuum method 
and in many cases is better, i.e. +0.001 per cent. absolute for low concentrations 
of CO,, and +0.25 per cent. of the CO, present for high concentrations. 


An extension of the above method has been developed for the estimation 
of CO. The gas, after the removal of the CO,, is passed over CuO at 350°C., then 
through another trap immersed in liquid air and finally, to the vacuum pump 
through another capillary restriction to contro] the flow. The CO, produced by 
oxidation of the CO is measured in the same way as the original CO,. This rapid 
method could be developed further to estimate other carbon-containing gases, but 
it still suffers from the disadvantage that hydrogen cannot be measured. 


3.3. pH-meter method 

This method, although developed before the rapid vacuum method was available 
for CO, determination, has now been largely replaced by it. Briefly, the principle 
of the method is to bubble the exhaust gas through N/1000 NaHCO, solution in 
which is immersed a glass electrode. This solution forms one arm of an electro- 
lytic cell, the other arm being a standard calomel electrode. The pH value of the 
NaHCO, solution, which is a function of the partial pressure of the CO, in the 
exhaust gas, is determined on a pH-meter. A calibration curve of pH value against 
CO, partial pressure can be obtained by passing known amounts of CO, 
through the cell. 


The method has the advantages of being rapid and continuous in operation. 
The main disadvantage is that the calibration occasionally changes because of 
electrical instability inherent in the present design of the cell system. 


The accuracy of the method is +2 per cent. of the actual amount of CO, in 
the sample under test, i.e. for a sample containing 2.0 per cent. CO, the absolute 
accuracy would be +0.04 per cent. The instrument is thus capable of determining 
small amounts of CO, (i.e. in combustion loss) to a high degree of absolute 
accuracy but it is not as suitable for larger amounts. 


For determination of CO, combustion efficiency the time required is 5 to 10 
minutes. A more detailed analysis can be made (i.e. for CO, CH,, CO,, etc.) by 
having three or more separate solutions and converting in turn each of these gases 
to CO,, as in the vacuum method. 


3.4. Other methods 


Other methods are either being developed or are under consideration, and the 
following notes are intended to indicate the potentialities of each method. 
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(i) The Poole Combustimeter 

This instrument is essentially a Wheatstone bridge, with a catalytic wire 
as resistance in one arm and a compensating non-catalytic wire as resistance in 
another arm. Provided that the temperature of the catalytic wire exceeds 
600°C. all combustible gases will be oxidised at the surface, but if the 
temperature is below 600°C., the method may be made selective for various 
gases. The heat of combustion released at the wire throws the bridge out of 
balance by an amount proportional to the concentration of “combustible” and 
its heat of combustion. The main interest in the instrument is that it furnishes 
an accurate method for the estimation of hydrogen. Provided that the 
temperature of the filament is kept below 300°C. the only constituents of 
combustion chamber exhaust gases which react are H, and CO and, up to 
the present, it has been impossible to separate the response of these two 
gases. Nevertheless, as CO can be estimated by an independent method the 
hydrogen can be found by difference. This will not affect the accuracy for H, 
as the instrument is much less sensitive to CO. The instrument is capable of 
giving a full-scale deflection for 0.88 per cent. hydrogen. This corresponds to a 
probable accuracy of +0.005 per cent. assuming the scale can be read to half 
a division. 


The great advantage of the method is that it is rapid and continuous. One 
disadvantage is that the apparatus must be initially calibrated, which means that 
provision must be made for preparing mixtures of hydrogen and air with 
high accuracy. 


(ii) Conductivity methods 
Some work has been done on the development of a conductivity method 
for the estimation of hydrogen but it was found that for the accuracy required 
(+0.005 per cent.) the instrument zero was insufficiently stable. This was 
probably because of the difficulty of maintaining the temperature of the 
incoming gas flow within a few hundredths of a degree centigrade. 


As the combustimeter is arranged to overcome this, it has now superseded 
the conductivity method in the present development programme. 


(iii) Infra-red gas analyser 

An instrument which appears to have great possibilities in this field is the 
infra-red gas analyser. In this instrument infra-red radiation is passed through 
a cell containing the sample under test and then to a detector cell containing 
only the constituent which is to be determined. A second identical beam of 
infra-red radiation passes through a second system, identical to the first except 
that the cell which contains the sample is replaced by one containing a standard 
non-absorbing gas (air free of CO,). The amounts of radiation received by 
the two detector cells will obviously differ by an amount which is equal to that 
absorbed by the relevant constituent in the sample and therefore proportional 
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to the concentration of this constituent. This differential effect is amplified and 
measured by an appropriate electrical circuit. The advantages of this instrument 
are its great accuracy (full-scale deflection for 0.03 per cent. CO, if necessary), 
its almost instantaneous response and its continuity. The gases in the exhaust 
which may be determined by means of this instrument are CO,, CO, CH, and 
other hydrocarbons. Hydrogen, not absorbing in the infra-red region of the 
spectrum, is of course not estimated. 


Making allowance for the full development of these various methods, it appears 
probable that analysis of the exhaust gas in determinations of combustion efficiency 
will be best effected by the use of a combustimeter in conjunction with the infra-red 
gas analyser. This will enable an analysis to be performed accurately, continuously, 
and rapidly. 


In this review, mention should be made of one other method which has been 
used extensively, namely the National Gas Turbine Establishment gravimetric 
method. This method uses fractional combustion at various controlled temperatures 
to convert each constituent to carbon dioxide and water. These products are then 
estimated gravimetrically by absorbing the CO, with soda asbestos and the water 
with anhydrone (magnesium perchlorate) and measuring the increase in weight. 
The various stages of combustion are similar to those used in the Cambridge vacuum 
method, the essential difference between the two techniques being in the method 
used to measure the combustion products. 


Further: details of the gravimetric method may be obtained from the 
reference quoted’). 


4. SAMPLING AND INTERPRETATION OF RESULTS 


With adequate analytical methods available, the most important problem 
remaining lies in extracting a sample which shall be representative of the exhaust 
gases; unless this requirement is satisfied, the overall accuracy of combustion 
efficiency assessment will suffer. 


In assessing efficiency it is the heat release within the combustion chamber 
which must be measured; any delayed energy release (the simplest case of this being 
a long flame), occurring beyond the outlet section of the chamber normally coinciding 
with the turbine blades, does not contribute to the most efficient working of the 
engine, and hence should not be included when computing efficiency. It follows 
that, in sampling, any delayed reaction occurring in the combustion chamber efflux 
gases should be arrested at the appropriate point; this can be accomplished for all 
practical purposes by chilling, or reducing gas temperature to a value where reaction 
is impossible, or at least extremely slow. 


Two courses are open, then, in obtaining the sample. Either some means must 
be provided for withdrawing at the outlet of the chamber, and immediately chilling, 
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a number of samples which when mixed provide a mean sample of the total exhaust 
gases; or the whole of the exhaust gases must be chilled on leaving the combustion 
chamber and sampled at some point downstream where they have mixed and are 
uniform in composition. Both practices have been followed at various establishments 
and each has its own particular merit and disadvantage. 


Of equipment for sampling directly at the combustion chamber outlet, perhaps 
the most successful design in current use is that developed at the National Gas 
Turbine Establishment’). This consists of a water-cooled sampling manifold 
pierced with a number of holes through which the sample is metered, the whole 
manifold being swept at constant angular velocity back and forth across the outlet 
section of the combustion chamber. By suitable adjustment of the sampling 
conditions the amount of gas admitted through any hole is made proportional to the 
local mass flow of gas, and thus the device supplies a sample truly representative of 
the areas swept by the holes. This equipment has proved satisfactory in operation 
where concentration gradients across the section sampled are relatively small, but 
where such is not the case, the fact that in effect the sample is drawn from a finite 
number of narrow arcs limits its performance. A logical development of this 
apparatus might be a single sampling tube which is made to scan completely the 
combustion chamber outlet (cf. television practice); mechanical difficulties may be 
discouraging but the construction of such equipment would be amply justified in 
some circumstances. 


Where gas composition varies widely across the outlet of the chamber perhaps 
the only satisfactory solution of the sampling difficulty lies in thoroughly mixing 
at first the whole of the exhaust gases and then withdrawing suitable quantities for 
analysis from a single sampling tube mounted at a convenient point downstream. 
As already stated, under certain conditions the use of this technique necessitates the 
chilling of the whole of the exhaust gases at the outlet of the chamber on test, and 
this problem has been tackled by injecting water from a grid of perforated pipes 
directly into the gas stream at a plane normally occupied by the turbine blading. 
The main problem involved is the design of a suitable gas mixer; the type in use at 
Thornton Research Centre was the outcome of trial and error; duct sections, baffles, 
etc., were modified and changed until gas composition across the exhaust pipe was 
found to be substantially uniform. A diagram of the sampling = finally adopted 
is given in Fig. 3. 


Both the methods of sampling outlined above suffer from one important 
limitation—where unburnt liquid fuel is present in the exhaust a truly representative 
sample cannot be collected. Thus the procedure followed for the practical 
determination of combustion efficiency varies according to the nature of the exhaust, 
i.e. whether it is entirely gaseous, or mainly gaseous but with some liquid fuel present. 
Experience has shown that at efficiencies below 85 per cent. the combustion loss 
will often include liquid fuel, and in order that the appropriate sampling procedure 
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THERMO-COUPLE SAMPLING TUBE 
spray THROTTLE VALVE r MIXER 
U 

COMBUSTION CHAMBER CIRCULAR BAFFLES 

CONCENTRIC 
OUCTING 
Fig. 3. 


Gas sampling gear at Thornton Research Centre. 


may be selected a rough approximation of efficiency is required. This can be 
obtained by estimating the temperature rise in the combustion chamber by means of 
a single thermo-couple mounted in the exhaust of the gas mixer (see Fig. 3). In 
view of the distinct differences in the techniques applied for the two exhaust 
conditions, it is proposed to deal separately with each. 


4.1. The determination of efficiency when exhaust products are 
wholly gaseous at the sampling temperature 

It is under these conditions that gas analysis methods of efficiency measurement 

are most accurate. Errors are governed only by the precision of analytical methods, 

the extent to which they are used, and the care taken in ensuring that the sample 

is truly representative of the exhaust gases; each of these factors, to a large extent, 

can be adjusted to suit the requirements of the test in hand. No measurements of 

combustion chamber air or fuel flows are necessary for the calculation of efficiency, 
because all data required are obtained from the analytical results. 


As already stated, the greater part of the combustion loss in conventional 
aircraft-type chambers will be accounted for by CO and unburnt hydrocarbons and, 
by neglecting the loss from any free hydrogen present in the exhaust (a step often, 
but not always, justified), two simple approximations to combustion efficiency may 
be determined; these are CO, combustion efficiency, already defined, and a closer 
approximation to true efficiency in which allowance is made for the liberation of 
heat in the formation of CO. 


CO, combustion efficiency may be determined by taking two samples, one (a) 
of the untreated exhaust gases, and one (b) of the exhaust gases treated in such 
a manner as first to remove CO, resulting from combustion in the test chamber, and 
then to oxidise unburnt carbon-containing constituents to CO,. It is usually 
convenient to install apparatus for treating sample (b) adjacent to the sampling gear, 
and a diagram of the necessary equipment is given in Fig. 4 with explanatory notes. 
Both samples (a) and (b) are examined for CO, content by one of the methods 
already outlined, the choice being governed by the speed and accuracy required. 


Then if CO,’=CO, in sample (a) expressed as percentage by volume of dry gas 
and CO,” =CO, in sample (b) expressed as percentage by volume of dry gas 
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CO, combustion efficiency = CO,’ + C0,” | x 100 per cent. . (1) 
= | 1- CO,’ CO,” | x 100 percent. . (2) 
In (2) ae ~—, x 100 per cent—the CO, combustion loss—may have 


resulted from either CO, hydrocarbons, or both, and hence cannot be interpreted 
directly as a thermal or energy loss. The closer approximation to true combustion 
efficiency, allowing for the proportions of CO and hydrocarbons in the combustion 
loss, can be determined by following the same sampling procedure as for CO, 
combustion efficiency but by analysing sample (a) for CO, as well as for CO, 
content, again by one of the methods described. 

Then if CO,””’=CO, produced by oxidation of CO in sample (a), expressed as a 

percentage by volume of dry gas: — 


(CO,” —CO,”’)+ k,CO,” 


i i pent. 3 
Combustion efficiency E CO,’+ C0,” ] x 100 percent. (3) 


Where (CO,” —CO,””) is proportional to the loss from hydrocarbons, and k,CO,”” 
is proportional to the loss from carbon monoxide. 


The combustion loss ! CO,’+ CO,” |. may be interpreted as an 


energy loss. 


k, is a constant for any particular fuel and is calculated as shown in Appendix I. 


The assumptions made in representing combustion efficiency by the expression 
(3) are that combustion loss from free hydrogen in the exhaust is negligible, and that 
the mean calorific value of the unburnt hydrocarbons (which may include cracked 
fuel, fuel residues, methane, etc.) is equal to the calorific value of the fuel burnt. 


A practical comparison between CO, combustion efficiency and the closer 
approximation to true efficiency is afforded by Table I, which shows typical results 
of tests on an aircraft gas turbine combustion chamber. 


For the most accurate measurement of efficiency consistent with reasonable 
rapidity, and neglecting only the small losses caused by intermediate combustion 
products and free carbon where present, hydrogen must be added to the exhaust 
constituents estimated. This can be accomplished by the use of one of the rather 
slow classical methods of gas analysis (e.g. Haldane) or by the gravimetric technique 
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used successfully at the National Gas Turbine Establishment). Rapid methods 
are under development at Thornton and particular attention is being given to the use 
of the Poole combustimeter for this purpose. 


Where, then, H, is the hydrogen content of the exhaust, expressed as a per- 
centage by volume of the dry exhaust gases : — 


_ k,H,.+k,CO,” +(CO,” —-CO,””) 
CO,” +CO,’ 


For the calculation of k, and k,, constants for any particular fuel, the reader 
is referred to Appendix I. 


Combustion efficiency = [1 ] x 100 per cent. 


4.2. Determination of combustion efficiency when combustion 
chamber exhaust products contain liquid fuel 
Under these conditions it is virtually impossible to obtain a representative sample 
of the exhaust because of the non-uniform distribution of liquid fuel—some may be 
in the form of widely-spaced droplets in the exhaust stream and some may be 
dribbling along the walls of the combustion chamber. Thus any analytical results 
from sample (b) would be meaningless. Fortunately, the presence of liquid fuel will 
not affect the sampling of the gaseous constituents of the exhaust, and CO,, CO 
and H, may still be measured accurately. The air/fuel ratio at which the chamber 
is Operating, or the total potential CO, in the exhaust (CO’,+CO,”), can no longer 
be calculated directly from the gas analysis data, and it becomes necessary to 
determine this ratio by the conventional method of measuring air and fuel flow 
rates separately. Every precaution should be taken to ensure accuracy in these 
measurements in order to avoid errors in the assessment of efficiency. 


To determine CO, combustion efficiency under these conditions the untreated 
exhaust gas sample is taken and examined for CO, content. 


Then if CO," is the theoretical percentage CO, content of the dry exhaust gases 
for complete combustion, calculated from the combustion chamber operating air / fuel 
ratio and the carbon/ hydrogen ratio of the fuel : — 

CO,’ 
CO, combustion efficiency = Co; * 100 per cent. 
2 

In determining the closer approximation to true combustion efficiency, taking 
into account the proportions of CO and hydrocarbons contributing to the combustion 
loss, sample (a) is taken and analysed for both CO and CO, as before. Whence :— 


(CO,'- CO,’ k,CO,”” 
CO,‘ 
Again, for maximum accuracy allowance must be made for any free hydrogen in 


the exhaust gases. Methods similar to those already suggested for liquid-free exhaust 
gases may be applied to a sample of untreated exhaust gas. Whence:— 


Combustion efficiency = [1 - ] x 100 per cent. 
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_k,H, + k,CO,’”+ (CO,'-CO,’-CO,””) 
CO,' 


Difficulty will be experienced sometimes in deciding whether or not the exhaust 
gases are free of liquid fuel. In such cases it is better to assume that they are and 
to sample accordingly; results of analyses will show whether sample (b) was valid, 
for if so CO,’+CO,” will be equal to CO,‘ within the limits of accuracy of the 
measurements taken. 

It has already been stated that the choice of analytical methods for estimating 
the important constituents of the exhaust will be governed mainly by the required 
speed and accuracy of efficiency determination. For guidance in this choice Table II 
has been prepared, which summarises limits of accuracy of the various methods 
considered, together with the time required for each determination. 


TABLE Il 


COMPARISON OF ANALYTICAL METHODS WITH RESPECT TO 
ACCURACY AND TIME REQUIRED 


Combustion efficiency =| 1 100 per cent. 


Accuracy and time required for determination of: 


| 
| co | Co, | H, | 
| 2 per cent. of | 2 per cent. of | 
pH | recorded CO | recorded CO, | 
| 5-10 minutes | 5-10 minutes 
| 
_— | At least as good as +0.005 per cent. | 
ae d) | absolute. Continuous 
|—- 
| +0.005 per cent. 
Combustimeter | | absolute 
| | Continuous 
Cambridge | +0.005 per cent. | +0.005 per cent. 
| 40 minutes | 50 minutes | 
Modified | +0.005 per cent. | +0.005 per cent. | 
absolute absolute | 
| 15 minutes* 20 minutes* | 


* N.B.—By using two sets of measuring apparatus, the CO and COs can be measured 
simultaneously by the modified method. The time for an analysis for CO and COs is 
then 25 minutes. 


183 


L. J. RICHARDS AND J. C. STREET 


5. POSSIBLE FUTURE APPLICATIONS 

The techniques described have been devised primarily for the purposes of 
laboratory rig tests where such equipment as gas chilling gear, mixing ducts, 
traversing samplers, etc., can be readily accommodated. Nevertheless, the methods 
have some possible applications in full-scale engine and field experiments where, 
in any case, other means of assessing combustion efficiency are often impracticable. 
In many instances, chilling of the exhaust gases at the plane of the turbine blades 
will be unnecessary, for in a well-developed combustion chamber chemical reaction 
beyond the turbine will rarely be observed; further, the passage of the gases through 
nozzle and turbine blading and thence into the jet pipe results in such turbulence that 
the prospects of producing a reasonably uniform gas mixture in jet or exhaust pipe 
are good. Thus, on engine test beds equipped with air and fuel measuring apparatus, 
methods similar to those already described could be adopted for efficiency 
determination, samples being taken from a number of stations in the engine exhaust 
pipe and subsequently mixed. In the field, and where exhaust products are wholly 
gaseous, the apparatus required for collecting samples necessary for efficiency 
assessment will be relatively simple and easily fitted to experimental aircraft. 
Unfortunately, under flight conditions giving rise to very low combustion efficiency 
and large quantities of unburnt fuel in the exhaust, sampling difficulties at present 
complicate the application of gas analysis as a means of efficiency determination, 
for under such conditions fuel and air flow rates cannot be measured readily 
with accuracy. 
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APPENDIX I 


Calculation of combustion efficiency from exhaust gas analysis 


Let percentage by volume of CO in exhaust gases =a 
Let percentage by volume of CO, in exhaust gases = 
Let percentage by volume of unburnt hydrocarbons in 
exhaust gases (expressed as CO.) =o 
Let percentage by volume of H, in exhaust gases =@ 


Consider a volume of exhaust gas equal to 22.4 litres. 


Volume of CO (at normal temperature and pressure) = 100 x 22.4 litres 
_ 28a 

Weight of CO = 700 stammes 

Similarly, weight of CO, = sed grammes 
44c 

Weight of CO, from hydrocarbons = 7009 Stammes 

Weight of H, = < grammes 

Therefore total weight of carbon = = erate) grammes 

and total weight of fuel = Ks Crete) (1 +A) grammes, 


where / is H/C ratio of the fuel by weight. 


Similarly weight of unburnt fuel = Ve (1 +A) grammes 
(assuming unburnt hydrocarbons are unburnt fuel). 


-. If Ky, Ky, Koo are the heats of combustion of the fuel, hydrogen and carbon 
monoxide, respectively, expressed as calories per gramme (or CHU/Ib.), 


2 2 
20 Kat Koo + (1+h) Ky 
Combustion loss = 12 x 100 per cent. 
100 (a+b+c)(1+h) Ky 


12 (a+b+c)(1+hA) Ky 


x per cent. 


Combustion efficiency = 100 — (combustion loss). 


Various approximations to the above expression can be used, as mentioned 
in the text. 
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(a) Neglecting hydrogen 


Combustion loss becomes : — 


28aK oo + 12c (1+h) Ky 
ate 
_ 120 
= x 100 per cent. 


In the method of sampling used in this report, 
If CO,’ ‘= percentage by volume of CO, in the untreated sample 
CO,” =percentage by volume of CO, in the treated sample 
CO,”” = percentage by volume of CO, from CO 


CO,‘ =total theoretical percentage of CO, for complete combustion 
derived from air/fuel ratio. 


Then when exhaust products are gaseous at sampling temperature : — 


_ kCO,” +(CO,” —CO,”) 
Combustion loss = (CO,”+C0,) x 100 per cent. 


_ __28Koo 
where k= 


When exhaust products contain liquid fuel: — 


kCO,”” 


Combustion loss = 


x 100 per cent. 


(b) Neglecting hydrogen and assuming CO is a total loss (i.e. CO, combustion 
efficiency) 


When exhaust products are gaseous : — 


Co,” 


CO, combustion loss = CO,” +00,’ x 100 per cent. 
In the presence of liquid fuel: — 
CO, combustion loss = aoe x 100 per cent. 
2 
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A SIMPLIFIED THEORY OF 
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“SIMPLE WAVES” 


A. GHAFFARI, Ph.D. 
(The University, Teheran) 


SUMMARY 


This note is concerned with a general and exact solution of the two-dimensional 
equations of motion of a compressible fluid, which is steady, irrotational and isen- 
tropic. The solution arrived at is the generalisation of the Prandtl-Meyer celebrated 
“wedge solution” to include an expansion around any smooth (or abrupt) 
convex boundary. 


We shali obtain the solution of equations for the most general type of motion in 
which the Cartesian components of velocity are connected by some functional 
relation throughout the field of flow. 


Taking into account the conditions of continuity, zero rotation and Bernoulli’s 
equation, together with the fundamental assumption, which implies that the 
components of fluid velocity are each a function of fluid speed only, it is found that 
the motion is necessarily supersonic. 


It is shown that the lines of equal speed, pressure (isobars) and density are all 
straight lines, although they are not necessarily concurrent as in the original Prandtl- 
Meyer theory, and further that the component of fluid velocity perpendicular to the 
isobar is equal to the local velocity of sound. 


The extended Prandtl-Meyer solution for the flow past an arbitrary fixed 
boundary is obtained analytically. 


This note is concerned with the plane motion of a compressible fluid, which is 
steady, irrotational and isentropic. It is proposed to obtain the most general type 
of flow in which there is a functional relation connecting the Cartesian components 
of velocity u and v, valid throughout the field of flow. 


This type of motion has been discussed in papers by Lees’) and Bateman”, 


[The Aeronautical’) Quarterly, Vol. 1, August 1949) 


Paper received October 1948. 
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Fig. 1. 
= +¢-—@ is the Mach angle, AB is the stream line, PQ is the isobar, PC is perpendicular to the 


isobar, w is the fluid velocity. 


but neither investigation brings out the relation of the flow to the well-known 
“wedge-solution” of Prandtl and Meyer“. 


To facilitate comparison with Meyer’s solution (as represented in Ref. 3), 


the notation of Fig. 1 is used so that 


u=—wsin< | 
v=wecos¢, | ( 


where w is the fluid speed and ¢ is the angular deviation of the flow from a direction 
parallel to the y-axis, measured in a clockwise direction. 
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A SIMPLIFIED THEORY OF “SIMPLE WAVES” 


The assumed functional relation between u and v is equivalent to the assumption 
that the velocity components u and v are each functions of the single variable w. 


NOTATION AND FUNDAMENTAL EQUATIONS 


Notation 


x,y rectangular co-ordinates. 
6 slope of isobar. 
9 angle between resultant fluid velocity and y-axis, 
measured in a clockwise direction. 
x 
p density of fluid. 
p static pressure. 


F,G functions defined for equations (10), (11) and (12). 


u,v components of fluid velocity parallel to rectangular axes of x and y. 

w resultant fluid speed = (u* + v*)?. 

? potential function. 

v stream function. 

y adiabatic index = 1.4 (for air). 

a local speed of sound. 

a, speed of sound at stagnation point. 

Wi, maximum which can be attained by fluid speed when p=0, p=0, 
and a=0. 


Fundamental equations 


The fundamental equations are given in the form 
(pu),+(pv),=90 (equation of continuity), . (2) 
u,=V, (condition of no vorticity), . (3) 


using suffixes, here and elsewhere, to denote partial (or total) derivatives. 


The standard pressure-density relation for a perfect gas is not assumed until 
the end of this investigation (i.e. from equation (13)). The main results are 
independent of this special assumption. 
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Since Bernoulli’s equation can be written in the form 


a’ dp —wdw . ‘ (4) 
p 
the equation of continuity (2) becomes 
UrtVy (up, + Vpy)/p=(uw, + Vw,) w/a’. 
whence + = (uw,.+ vw,) (w/a’), 
or W, (Uy — uw /a’)+Ww, (Vy — vw/a?)=0. ; (5) 


The condition for irrotational flow (3) becomes 
Hence, from the two equations (5) and (6) 


Uw (Uy — UW/ a?) + Ve (Ve — vw/a*)=0, 


or 
Ur? + Vip? = (Ul + (w/a?) 
[d (u? + v?)/dw] (w/a?)=w?/a’. . (7) 
Now from (1) 
Uy, = — SIN ¢ —W COS ¢ 
(8) 
and Vy = COS ¢ —W Sy SIN ¢ 
Therefore 


Uy? + Ve? = 1+ we’. 
Hence from (7) 

=w?/a’ 
or 


dw 


y =(w?—a’)/w’a?, ‘ 


It is clear from the last equation that w? > a’, that is to say w cannot be less 
than a, i.e. the motion is supersonic everywhere. 
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Before integrating equation (9) it is convenient to determine the 
features of the motion. If we write 


F o+xu+yv, 


where @ is the velocity potential, then 


brt+Ut+XU, + W,, 
and Fy = —by t+ V+XUy + =(XUy + Wy. 
0 (F, w) 
Therefore F,w,—F,w,=0, or 0, 
d(x, y 


i.e. F is a function of w only, say F (w). 


dF (w) 
Hence - W, U+(XUy + We, 
dw 
dF (w) 
or — =XUy + Ww, 
dw 
or C=Ax+By, 


where A, B and C are functions of w only. 


main 


(10) 


If w =constant, the equation (10) represents a straight line. The curves of equal 
fluid speed w are also the curves of equal density and pressure, i.e. they are isobars. 


The isobars therefore form a system of straight lines. 


If the slope of an isobar is denoted by @, then 


sin 4 au,,/w, from (7) 


and cos 6=av,,/w, 


the signs being chosen so that @ is positive when ¢ is small and positive and increases 
with w. Hence the component of the fluid velocity perpendicular to the isobar is 


w cos (¢—<)=vcos 6—u sin 6=(vv, +uu,)a/w=a, 


i.e. it is equal to the local speed of sound (Fig. 1). 
Similarly, if we write 
G= —v+xpv—ypu, 
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Vet pV +x (pu), 
(pV )w — (PU)w] 


G,= —vy— put+x (pv), —y (pu), 


=[x (pv). y We. 


0(G, w) 
Therefore G.w,—G,w,=0, or 


i.e. G is a function of w only, say G (w), 
and v= xpv—ypu—G (w), 
or v= pw (x cos ¢+y sin (w), 
The equation of the stream lines is therefore 
x cos ¢+y sin ¢=R 
where R (¢) is some function of ¢. 
R (¢)= [G (w)+ constant] / pw. 


Hence, if any one stream line is known, as for example a solid boundary of the 
fluid, then x, y, ¢ are also determined everywhere as functions of w. 


To avoid the ambiguity in sign which is implicit in (9), we determine ;¢,, 
from the equation 


tan ¢= —u/v. 
Then 
= — (VU, — UV,,)/(v? sec? ¢)=(w/a)(v sin cos 4) (1/w*) 


(sin cos ¢ —cos sin 7)/a=sin 
We can complete the integration in the case of a perfect gas, for which 
a’ =a,’ —4(y—1)w’. 
We now express w and a in values of the angle 
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We have 


a=weosy, from (11) 


and, taking into account equation (13), we obtain 


(y— 1)+cos? x’ 
4 dw 2 cos yx sin x dy 
4 (y—1)+ cos? 


sin yd lw sin*\d 
\GX tan ¢ \aX 
a w (y— 1)+cos* x 


=-dy+ Ddx 


Therefore pe (A tan 


where A y 


since g and x vanish together. 
All the variables can now be expressed in terms of the angle y. For 
¢=(1/A) tan-' (A tan 
6=x+¢, 
tan y= (1/A) tan (A), 


2a,” 


= 4w,,? [(1 +A?)— (1 —A?) cos 2A8], 
u=—wsin¢, 


a=W COS x. 


These expressions for the components of velocity will be recognised to be 
identical with those of the Prandtl-Meyer “wedge-solution” with one significant 
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modification. The angle 6 has not been assumed to be the ordinary radial 
co-ordinate. Hence the isobars (6=constant) are not necessarily concurrent straight 
lines as in the original Prandtl-Meyer solution. 


We have in fact obtained analytically the extended Prandtl-Meyer solution for 
the flow past an arbitrary fixed boundary. To establish this we have only to show 
that one stream line can be prescribed arbitrarily in our solution. 


Now let a stream line be given in the form 
x=X(%), y=y(>) 


¢ being the direction of the tangent to the stream line. 


Then the equations (14) and (16) determine w, a and therefore p, in terms of ¢ 
along the stream line. Equation (12) then determines the unknown function G (w), 
and ail the stream lines are determined. 


The preceding analysis therefore furnishes a formal proof of the extended 
Prandtl-Meyer solution, which has hitherto been adopted as an intuitively obvious 
generalisation of the original “wedge-solution.” 


The author wishes to thank Professor G. Temple, F.R.S., for his 
valuable suggestions. 
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